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Аннотация 

Введение. Задача анализа устойчивости пластин и оболочек в условиях ползучести актуальна для элементов 
конструкций из материалов, обладающих свойством старения, находящихся под действием длительных нагрузок, 
поскольку потеря устойчивости может происходить резко и задолго до исчерпания прочностного ресурса 
материала. Вопросы совместного учета геометрической нелинейности и ползучести в задачах выпучивания 
пластин в настоящее время остаются слабо изученными, существующие программные комплексы не позволяют 
выполнить такой расчёт. Целью настоящей работы выступает разработка алгоритма расчета на устойчивость 
прямоугольных пластинок с начальной погибью, испытывающих действие нагрузок в срединной плоскости с 
учетом геометрической нелинейности и ползучести. 

Материалы и методы. При получении разрешающих уравнений в основу положены геометрические и 
статические уравнения теории гибких упругих пластин. Физические уравнения выводятся из предположения, что 
полные деформации равны сумме упругих деформаций и деформаций ползучести. Окончательно задача была 
сведена к системе из двух дифференциальных уравнений, в которых в качестве искомых функций выступают 
функция напряжений и прогиба. Решение полученной системы уравнений выполнялось численно с помощью 
метода конечных разностей в сочетании с методом последовательных приближений и методом Эйлера. В 
качестве граничных условий для функции напряжений используется рамная аналогия, как в случае плоской 
задачи теории упругости. 

Результаты исследования. В рамках поставленной цели разработан алгоритм расчета и представлено решение 
задачи для пластины, сжимаемой в одном направлении равномерно распределенной нагрузкой. Исследован 
характер роста перемещений при различной величине нагрузки и начальной погиби. Установлено, что при 
достижении вертикальными перемещениями величин, соизмеримых с толщиной пластинки, скорость их роста 
начинает затухать даже при нагрузке больше длительной критической. 

Обсуждение и заключение. Результаты анализа устойчивости с использованием разработанного алгоритма 
показывают, что рост прогиба пластины при рассмотренных граничных условиях ограничен, потеря 
устойчивости не наблюдается при любых значениях нагрузки, не превосходящих мгновенную критическую. Это 
говорит о возможности длительной безопасной эксплуатации таких конструкций при нагрузке менее мгновенной 
критической. 


Ключевые слова: устойчивость, ползучесть, пластина, геометрическая нелинейность, физическая нелинейность, 
начальные несовершенства, метод конечных разностей 
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Введение. Анализу устойчивости тонкостенных конструкций в виде пластин и оболочек уделяется большое 
внимание, поскольку такие конструкции широко используются в строительстве и других областях 
техники [1-3]. Одной из актуальных задач в области расчета пластин и оболочек является анализ их напряженно- 
деформированного состояния в условиях ползучести, что подтверждается значительным числом работ, 
опубликованных в последнее время по данной проблеме в отечественных и зарубежных источниках. 
Так, в работах [4-8] исследуются вопросы выпучивания при ползучести композитных тонкостенных 
конструкций. В статье [9] рассматривается задача устойчивости функционально-градиентных пластин с учетом 
зависимости свойств материала от температуры. В работе [10] в задаче выпучивания композитных пластин 
применены методы стохастического анализа. В статьях [11-17] изложены вопросы устойчивости вязкоупругих 
пластин и оболочек при воздействии динамической и следящей нагрузки, а в статье [18] рассматриваются 
пластины средней толщины с учетом зависимости свойств материала от времени. Возникающие при решении 
указанных задач математические трудности приводят к тому, что многие исследователи ограничиваются 
линейными законами вязкоупругого деформирования или рассматривают случай установившейся ползучести. 
Большие возможности в задачах расчета пластин и оболочек с учетом ползучести открывает метод конечных 
элементов. Однако в современные расчетные комплексы, такие как АМЗУ$, АБадиз, ЛИРА и др., заложен 
ограниченный набор реологических моделей, применимых для конкретных материалов в фиксированном 
диапазоне напряжений и температур. Существует необходимость в альтернативных методах расчета, 
подходящих для произвольных законов вязкоупругого деформирования, в том числе нелинейных. 
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Целью настоящей работы явилось построение системы разрешающих уравнений для задачи выпучивания 
пластин прямоугольной формы с нелинейными вязкоупругими свойствами при действии усилий в срединной 
плоскости, с учетом больших перемещений, а также алгоритма ее решения. Отметим, что проблема устойчивости 
элементов конструкций с учетом ползучести не может быть решена, как проблема чистой устойчивости. Ее 
решение требует наличия возмущений в виде начальных неправильностей. Обычно начальные несовершенства 
задаются в виде начальной погиби или эксцентриситетов приложения нагрузок. 

Материалы и методы. Методику расчета рассмотрим на примере пластинки с начальной погибью у(х,у), имеющей 
шарнирное опирание по контуру и сжимаемой в направлении оси х распределенной нагрузкой р [кН/м] (рис. 1). 


х 


Рис. 1. Расчетная схема 


В рассматриваемом случае при наличии ползучести, если сравнивать с теорией упругих гибких пластин, 
отличие будет проявляться только в форме физических уравнений. Полные деформации могут быть 
представлены в виде суммы деформаций срединной плоскости (переходящей в поверхность) и изгибных 
деформаций, которые вызваны изменением показателей кривизны срединной поверхности: 


02 02 02 
Е 8, =50-—2 у» = -22 я (1) 
о 0х’ * д” дхду 
где &;и &, — полные линейные деформации; у,— полные угловые деформации; = и = — линейные 


деформации срединной поверхности; у — угловые деформации срединной поверхности. 
Для деформаций срединной поверхности может быть записано уравнение неразрывности деформаций [19]: 


д2ео ‚ 928 02° _[ 02 (м+м) , 92 (+) 0? (и+щ) |[02м, д? дм, 


(2) 
ду? 0х? одхду дхду 9х? ду? дхду 9х? ду? 
Для материалов с вязкоупругими свойствами полные деформации можно представить в виде: 
1 ыы г ОМУ) пн 
=, = © —=У0, ) +=,;=, = = (о, —уУб, ) +5; 7, = й и (3) 
где 5..5, — деформации ползучести; ЕЁ — модуль упругости; у — коэффициент Пуассона; 0», ду, у — 


величины компонентов напряжений по соответствующим направлениям. 
Выразив в (3) компоненты напряжений через деформации, запишем физические соотношения в обратной форме: 


с. = (. + УЕ - (& +е;)}; в Е (: + УЕ (в +=); = =] ) (4) 
Ут > ” т Та ме а 2(1+у)`® ” | 
Связь внутренних силовых факторов с напряжениями определяется интегральными соотношениями: 
1/2 й/2 1/2 #/2 1/2 1/2 
№, = | 0,4;М,= [6,4;5= | 1,42; М, = [| 6,24;М,= | ©,24;Н = | 1,24, (5) 
—й/2 —1/2 —1/2 —й/2 —й/2 —1/2 
где №, и М, — погонные продольные силы; 5 — погонные сдвигающие силы; М, и М, — погонные 
изгибающие моменты; Н — погонные крутящие моменты; й — толщина пластинки. 
Далее подставим (1) в (4), а также (4) в (5). В итоге получим: 
Ей Ей 
„= (оу )- №№, = (оу )-№;5 = 7-5”; 
1-2‘ ° > 1-\2\: | 2(1+у) 
02%  02м 02 2 02% 
М +У М*; М, =-Б +У м:;Н =-Б(1-у) Н*, (6) 
х 2 2 х У 2 2 3? 
ох ду ду дх дхду 
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Ейз Е #2 
где = — цилиндрическая жесткость пластинки, № = | (= +У5* ) 2, 
12(1-у2) 14 У2 42° . 
Е #2 Е #2 
М№* = 5 + УЕ” | 42, 5* = * 2, 
7" 1-у2 в. й я 2 +у) 1: 
Е #2 Е #12 Е 2 
М' = = + УЕ” | 242, М* = = +УЕ* | 248, Н* = |1" 202. 
“ ]-У2 й ы ;) 7 1-\2 И 7 я 50 +у) 


Величины №, №, 5*, М!, М;, Н* имеют размерность внутренних усилий и определяют вклад деформаций 
ползучести в перераспределение усилий. 
Статические уравнения теории гибких пластин имеют вид [19]: 
ОМ, а 05 _ (05 ОМ, _ 
; 


+ 0; 
0х ду 0х ду 


9?М 02Н 0?М, 92 (и+м) 2 (и+ж) 92 (и+ и) р 
+2 + +М№, +М, +25 + 
дх? дхду ду? | дх? ду? дхду 
Здесь 4 — нормальная нагрузка по поверхности пластины, которая в данной задаче равна нулю. 
Удовлетворить первым двум статическим уравнениям можно при помощи функции напряжений Эри, 


введенной по формулам: 
2 2 2 
о ро В 0°Ф_ (8) 
о 0 г 9 Охду 


После подстановки последних трех равенств из (6) в последнее статическое уравнение в (7) и с учетом (8), 
получим первое разрешающее уравнение: 
д2ф 0? (+ и) д2ф 0? (ж+ и) 2 д2ф д? (ж+ и) 

+ 
9х? ду? ду? дх? дхду —дхду 
02М* _0Н* 0?М; 
где 4“ = +2 + : 
9х? дхду ду? 


Чтобы получить второе разрешающее уравнение, необходимо выразить из (6) деформации срединной 


(9) 


РУ =а+а“+ 


поверхности: 
М. — УМ, + №" —уМ* 2 2 
20 =—^ - ы В м УМ, |; 
ЕЙ ЕЙ ду? 0х? | 
№, УМ, + № — УМ: 2 р 
&0 =—^ й 2 ЕЕ +2 к: УМ" |; (10) 
7 Ей ЕЙ \. Ох? ду? ? 
2(1+ 2(1+ 2 
И | 
Ей ЕЙ дхду 
Подставив (10) в (2), получим: 
2 2 
Фе д? (ити) д? _ 0? (учи) 0? (+) 9 бщ 1 2(1+)2° + 
Ей дхду дхду дх? ду? 0х? 0? Ей дхду (1) 
02мм" 02М№, | дм 92М№ 
+У + ь : 
9х? ду? ду? Ох? 


Таким образом, для рассматриваемой задачи получена система разрешающих уравнений из двух 
дифференциальных уравнений четвертого порядка (9) и (11). Уравнения (9) и (11) являются нелинейными. В 
полученных уравнениях величины Ф и \ представляют собой функции координат х, у, а также времени (. В 
явном виде в данных уравнениях время отсутствует, зависимость от времени закладывается в деформации 


ползучести &,, =", у.,, которые учитываются введением интегральных величин №, №, 5", М!:, М;,Н*. 


Для изображенной на рис. 1 расчетной схемы граничные условия записываются в виде: 


2ф 2ф 
ОЕ — =-р;5 = О =0;и=0;М, =0; 
“ ду? дхду ° 
@Ф 2Ф (12) 
приу=0,у=Ь: М, = = 0;5 =— = 0;ж=0; М, =0. 
- д? дхду 
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Уравнение (11) при малых перемещениях, в случае пластины из упругого материала, представляет 
бигармоническое уравнение, которое используется для решения плоской задачи теории упругости в 
напряжениях. Граничными условиями по функции напряжений для бигармонического уравнения может служить 
рамная аналогия. Контур пластины рассматривается как рама и функция напряжений на контуре будет равна в 
ней изгибающему моменту М, а ее производная по нормали к контуру — продольной силе М. Эпюры М и №Мв 
раме могут быть построены в одной из основных систем метода сил (ОСМО). Основная система, а также эпюры 
изгибающего момента и продольной силы в раме показаны на рис. 2. 


Рис. 2. ОСМС и эпюры изгибающих моментов и продольных сил 


Если вертикальные перемещения не превышают четверть толщины пластинки, то можно принять усилия в 
срединной поверхности не зависящими от координат хиу (М, =-р,М№, =5=0) и использовать для расчетов 


линеаризованное уравнение: 

02% 02% 

РУ\и+р =4"-р-—*. 
Ох? )х2 


(13) 


Аналитическое решение системы уравнений (9) и (11) связано с большими трудностями. Авторы предлагают 
решать данную систему численно. Используется метод конечных разностей (МКР) в комбинации с методом 
последовательных приближений. Для определения деформации ползучести во временной области применяется 
метод Эйлера. В качестве первого этапа выполняется решение для упругой пластинки. Нагрузка р 
прикладывается ступенчато с небольшим шагом. При начальных значениях нагрузки прогибы и! вычисляются 
путем решения упрощенного уравнения (13). Затем выполняется подстановка вычисленных значений у в 
дифференциальное уравнение (11). Это позволяет определить функцию напряжений. Следующим шагом 
является решение дифференциального уравнения (9) с использованием известных величин функции Ф, что 


позволяет определить узловые величины прогибов и”. После этого выполняется подстановка значений 
и, = (и + ) /2 в уравнение (11). Повторение итерационного процесса на каждом шаге происходит до тех пор, 
пока относительное расхождение с нормами векторов узловых значений прогибов и, и и” больше заданного 


значения (авторами оно принималось равным 0,1 %). Для второго шага по нагрузке начальным значением у" в 


каждом узле выступает конечный результат, полученный на первом шаге. Методика расчета во временной 
области с учетом ползучести аналогична. Интервал времени, на котором исследуется процесс, делится на шаги 
ДЕ. В случае задания закона вязкоупругого деформирования в дифференциальной форме величины деформаций 
ползучести на шаге #+ А! вычисляются на основе известной скорости их роста в момент времени Ё с 
использованием аппроксимации по Эйлеру: 


* 


0Е 
=, =8, +-—— ДЕ. 14 
ИА 1 6 ( ) 


Блок-схема алгоритма расчета на ползучесть приведена на рис. 3. 

Отметим, что система уравнений (9) и (11) позволяет использовать схемы более высокого порядка точности, 
например, метод Рунге-Кутта четвертого порядка. При этом для достижения той же точности результатов можно 
задавать заметно большие шаги по времени. Однако при увеличении шага есть вероятность не уловить эффекты 
неустановившейся ползучести в начальные моменты времени. А при одинаковом шаге по времени метод Рунге- 
Кутта, по сравнению с методом Эйлера, требует выполнить в четыре раза больше операций. 
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Начало 


Решение упругой задачи: 
а ь— 


Определение уз”, из ур-я (9) 
нет 


и, =(и и, ,)/2 


Подстановка з», в (11) и определение Ф, 


пах. 
—№ 


а 100% 


05* 
числение с, — 
0! 


ВЫ 
нет 


Вывод результатов 


Рис. 3. Блок-схема алгоритма расчета на ползучесть 


Результаты исследования. Рассмотрена полимерная пластинка из поливинилхлорида размерами а = 2 М, 


= Б=2 м, р =1см при Е = 1480 МПа, у = 0,3. В качестве закона, определяющего скорость роста деформаций 
Е ползучести, было принято нелинейное уравнение Максвелла-Гуревича: 
з Е; 9 а _ О 
- д 9 1 х, У, и х, У, 
В Е т 
5 3 
> м: *. 
= Л =5(°, -,5,)-Е.8}; (15) 
8 
= А 3 

т" — то с | и ра ( И бу ) Е„5,, 2 

т 2 тах 
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где 6, — символ Кронекера; о, = Л /З, Л =о, +в, — первый инвариант тензора напряжений; Пр, Е» ит — 
реологические параметры материала, называемые начальной релаксационной вязкостью, модулем 
высокоэластичности и модулем скорости. 

Индексами гг в формуле (15) обозначено направление главных напряжений. 

Для ПВХ значения реологических параметров авторами были взяты из работы [20]: Е»= 5,99-103 МПа, 
т*=12,6 МПа, (у = 5,44.107 МПа-с. Форма начальной погиби уо(х,у) принималась в соответствии с первой 


формой потери устойчивости пластинки из упругого материала: 
. ЛХ. ЛУ 
\ (ху) = ляп 5т- (16) 


Для пластины из упругого материала без начальных несовершенств критическая нагрузка, в случае целого 
отношения сторон а/б, определяется по формуле [19]: 
_ 412 
Р» = ГЫ 
Для верификации разработанного алгоритма расчета первым этапом было выполнено решение тестовой 
упругой задачи и сравнение результатов с расчетом в конечно-элементном пакете ЛИРА-САПР (рис. 4). 
Величина стрелы начальной погиби /о задавалась равной 0,15 мм. Размер сетки при использовании МКР 
составлял 20х20, количество шагов по нагрузке — 200. При расчете в ПК ЛИРА-САПР выполнялось разбиение 
пластинки треугольными конечными элементами с шагом триангуляции 0,1 м. Величина шага по нагрузке 
принималась такая же, как и при использовании МКР. Значение критической нагрузки для упругой пластинки, 
вычисленное по формуле (17), составило 1340 Н/м. В таблице 1 показано сравнение вертикальных перемещений 
в центре пластинки для различных значений нагрузки, полученных по авторской методике и при помощи метода 
конечных элементов (МКЭ). Прогибы, вычисленные с использованием двух альтернативных методов, 
достаточно близки, за исключением нагрузки 1330 Н/м. Отклонение при этой величине нагрузки можно 
объяснить тем, что при приближении к критической нагрузке перемещения устремляются в бесконечность. 


(17) 


ЕЕ ЕЕ Е ИИ ПОИИНИВ 
0 0,0365 0,0456 0,911 1,37 1,82 2,28 2,73 3,19 3,65 


Рис. 4. Изополя вертикальных перемещений в ПК ЛИРА-САПР ( р = 1330 Н/м): 
нелинейное загружение 1, изополя перемещений по 7(@), еденицы измерения — мм 


Таблица 1 


Сравнение результатов расчета по авторской методике и при помощи МКЭ 


и’-103, мм 
р, Н/м 
ЛИРА-САПР Авторская методика 

133 16 16 
266 37 37 
399 63 63 
532 98 99 
665 146 148 
798 218 221 
931 336 342 
1064 562 578 
1197 1163 1229 
1330 3646 4229 
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В статье рп показана возможность перехода от решения упругой задачи расчета пластин к решению в конце 
процесса ползучести. Величину длительной критической нагрузки р» можно получить, заменив цилиндрическую 
жесткость Р упругой пластинки на длительную цилиндрическую жесткость РБ», которая определяется по формуле: 

фе Е 
2 разв), (18) 
Е’ . И - * > 
3 - 

Для вязкоупругих стержней и круглых пластин было ранее установлено, что в случае р<р», рост 
перемещений во времени замедляется, и стрела прогиба приходит к конечному значению. Если р = Ро, прогибы 
растут с постоянной скоростью. При р > р» скорость роста прогибов возрастает. 

Авторами также проведен анализ характера роста прогибов во времени для р <р», р=р» и р> ро для 
различных значений максимальной начальной погиби /о. Кривые изменения прогиба с течением времени в центре 


где @ = 


пластины при р = 0,9 р»,р = р» ир=1,1 р» приведены соответственно на рис. 5-7. 


и, мм 
1,4 т т т т т т т т т 
——— 10,025 мм 
1,2 | — 1=0,05 мм 


—— /=0,1 мм 
1 —— ^=0,15 мм 


1,0 


0,8 


0,6 


0,4 


0,2 


0 2 4 6 8 10 12 14 16 18 20 час 


Рис. 5. Кривые изменения прогиба в центре пластины во времени при р = 0,9р» 


у, мм 


4,5 т т Т т т т т т т 
— — /№-=0,025 мм 


—— /=0,05 мм 
— /=0,1 мм 
—— /=0,15 мм 


4,0 
3,5 
3,0 
28 
2,0 
1,5 
1,0 


0,5 


0 2 4 6 8 10 12 14 16 18 20 час 


Рис. 6. Кривые изменения прогиба в центре пластины во времени при р = р» 
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Е 


у, мм 


8 т т т т т т т т т 


—— /=0,15 мм 
—— /=0,1 мм 
—— /—=0,05 мм 
= /=0,025 мм 


0 1 2 3 4 5 6 7 8 9 10 ББ час 
Рис. 7. Кривые изменения прогиба в центре пластины во времени при р = 1,1 р» 


Обсуждение и заключение. Из рис. 5 видно, что при р <р» стрела прогиба всегда приходит к конечному 
значению, независимо от значений начальных несовершенств. В то же время при р > р» полученный в [21] 
характер роста прогибов имеет место только при малых перемещениях. При достижении прогибами значений, 
превышающих примерно четверть толщины пластины, скорость роста деформаций начинает убывать даже в 
случае нагрузок, превышающих длительную критическую. Следует также отметить полное отсутствие участка с 
возрастающей скоростью роста перемещений для пластин с большими начальными искривлениями. Объяснить 
выявленные эффекты можно перераспределением усилий М№,, №,, 5 в срединной поверхности. 


Резюмируя вышесказанное, можно сделать вывод, что вертикальные перемещения пластины, шарнирно 
опертой по контуру, при действии сжимающей нагрузки по [611 (67 оси, всегда приходят к конечному значению, 
если нагрузка не превышает мгновенную критическую. Иначе говоря, при рассмотренном закреплении и 
нагружении пластина в условиях ползучести находится в устойчивом равновесии. 

Полученные уравнения и алгоритм вычислений позволяют рассчитывать пластинки из произвольных 
вязкоупругих материалов при любых вариантах закрепления. Закон связи между напряжениями и деформациями 
ползучести также может быть задан произвольно. 
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